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Abstract 



In the conventional QCD sum rule calculation for the ttNN coupling be- 
yond the soft-pion limit, there are three distinct Dirac structures, each of 
which can be used to construct separate sum rule. We point out subtleties 
in the previous sum rule results, based on one of the Dirac structure, asso- 
ciated with using either the PV or the PS coupling schemes in modeling the 
phenomenological side. We propose a sum rule coming from a different Dirac 
structure, which is independent of the coupling schemes used and has less 
uncertainty in the OPE. The obtained value is g n N = 9.76 ± 2.04 where the 
uncertainty mainly comes from the quark-gluon mixed condensate. 
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Since first introduced by Shifman, Vainshtein and Zakharov [||, QCD sum rule has been 
widely used to study the hadron properties such as masses or couplings of baryons @]. QCD 
sum rule is a framework which connects the physical parameters with QCD parameters. 
In this framework, a correlation function is introduced in terms of interpolating fields con- 
structed from quark and gluon fields. Then, the correlation function, on the one hand, 
is calculated by Wilson's operator product expansion (OPE) and, on the other hand, its 
phenomenological "ansatz" is constructed. A physical quantity of interest is extracted by 
matching the two descriptions in the deep Euclidean region (q 2 = — oo) making use of the 
dispersion relation. The extracted parameter therefore should be independent of the possible 
ansatz in order to be physically meaningful. 

One important quantity to be determined in hadron physics is the pion-nucleon coupling 
constant, g^N- Empirically its value is known to be around 13.4 but it is of interest to 
determine the coupling from QCD. QCD sum rule can be used for such purpose and indeed 
there are such calculations of g^N Reinders, Rubinstein and Yazaki [[J calculated 

g-irN by retaining only the first nonperturbative term in OPE. Later Shiomi and Hatsuda 
(SH) || improved the calculation by including higher order terms in OPE. SH considered 
the two-point correlation function for the nucleon interpolating field Jn, 

U(q,p) = % [ d 4 xe l ^(0\T[J N (x).I N (0)Mp)) , (1) 



and evaluated the OPE in the soft-pion limit (p M — > 0). 

More recently, Birse and Krippa (BK) |||| pointed out that the use of the soft-pion limit 
does not constitute an independent sum rule from the nucleon sum rule because in the limit 
the correlation function is just a chiral rotation of the nucleon correlation function, 

U(q) = i J rf 4 xe^(0|T[J^(x)J^(0)]|0) . (2) 

Therefore, BK considered the sum rule beyond the soft-pion limit. However, in their ap- 
proach, the OPE contains some condensates which are not precisely known. This causes 
large uncertainty in their result, g n N = 12 ± 5. 

In this letter, we construct an alternative sum rule for g n N which contains less uncertainty. 
We start by noting that the previous calculations have coupling scheme dependence, namely 
the dependence on either the pseudoscalar or the pseudovector type of interaction in the 
construction of the phenomenological side. We stress the importance of doing the sum rule 
calculation independent of the coupling scheme. To see the origin of the problem, consider 
saturating Eq.(|l]) with a nucleon intermediate state. Then we have, 

II(g,p) ~ \ N l±?L{N{q)\J N \*(p)). (3) 
qz _ m z 

The coupling, g nN , is defined to be the residue of the last term. 

— d — p* ~\~ 771 

(N(q)\J N \n(p)) ~ _ - m2 ^NgnN + non pole term (4) 

How to treat the non pole term will depend on the coupling scheme used and as we will see, 
influence the subsequent sum rule analysis for g nN 
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To be more specific, the pion-nucleon coupling has a unique definition as a matrix element 
of a pion field sandwiched by the nucleon states. This matrix element can appear in Eq. ([!]) 
by, with the use of LSZ reduction, moving the pion state into the time-ordered product and 
inserting on-shell nucleon intermediate state. This corresponds to the residue of the double 
pole while the rest contains the single nucleon pole which has dependence on the coupling 
scheme and therefore needs to be treated carefully. The calculation beyond the soft-pion 
limit is natural in this consideration. When starting from Eq. (lip, to leading order in the 
pion momentum, there are three distinct Dirac structures which can in principle be used to 
calculate g^N- We discuss them separately and propose a new sum rule to determine g n ^. 

Regarding the existing sum rule calculations of g n N @-f|, one question is whether the 
results have any dependence on the nNN coupling scheme adopted. This issue, within our 
knowledge, has not been properly considered but it can provide important aspect of how 
one can proceed in the construction of the sum rule. This question is not specific to 7rNN 
coupling and it can generally apply to any sum rule calculation involving the pion. 

In QCD sum rule, the phenomenological side of the correlator is usually constructed by 
making an "ansatz" based on some effective models. In the case of the two-point correlation 
function in Eq. (|l|), one often uses the pseudoscalar interaction Lagrangian, 

£ ps = g-KN^%lhT ■ 7T^ , (5) 

in constructing the phenomenological side. However, the pseudovector interaction La- 
grangian, 

Cpv = §^^m t ■ , (6) 

can also be used for constructing the phenomenological side. At present, we do not know 
which Lagrangian is more reliable in describing the pion interaction with nucleons. The two 
descriptions are equivalent when the participating nucleons are on-shell but they are usually 
not when the nucleons are off-shell. Since on-shell properties of a particle are extracted 
by matching the OPE with the phenomenological side in the deep Euclidean region, the 
pion coupling constant determined from QCD sum rules may have the coupling scheme 
dependence. Let us employ the two descriptions in the phenomenological side of the sum 
rule and see if they lead to any difference. 

Using the interaction Lagrangians in the two-point correlation function, Eq. (|l|), we find 
that there are four distinct Dirac structures, (1) 75, (2) 750^, (3) 75 p 1 , (4) 75 The fourth 
structure appears only in the pseudovector case with the residue factor, p 2 — 2p ■ q and this 
structure is definitely scheme- dependent. In principle, any of these structures can be used 
to calculate the coupling. Shiomi and Hatsuda (SH) [||] considered the 75 structure in the 
soft-pion limit. 

When the pseudoscalar interaction Lagrangian is used, the 75 structure without taking 
the soft-pion limit takes the form 

ig*N>?N H Wttn^nP ■ q , 7 v 

(q — p) 2 — m 2 + ie (q 2 — m 2 + ie) \{q — p) 2 — m 2 + ie] 

where is coupling of the nucleon interpolating field, J N , to the physical nucleon, m is the 
nucleon mass. The dots indicate the contribution from continuum whose explicit forms are 
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not necessary for our discussion. In the soft-pion limit, the double pole structure disappears, 
leaving only the simple pole which was used by SH. 

The corresponding expression in the pseudovector coupling scheme is 

ig nN X 2 N P 2 /2 + _ _ _ ^ 



(q 2 — m 2 + ie) [(q — p) 2 — m 2 + ie] 

A crucial difference from above is that now the correlator contains the double pole only 
More importantly, in the soft-pion limit, the correlator is zero and one can not extract any 
information about the coupling. Clearly, the method by SH is scheme-dependent. Another 
problem is that since SH considered the simple pole, the additional simple pole coming from 
N — ► N* transition || can not be separated from the nucleon simple pole. 

What will then be the reliable procedure using this 75 structure ? The common pole in 
both coupling schemes is the double pole. Consider expanding the denominator of the double 
pole in terms of the pion momentum and keeping only the leading term. The resulting term 
can also be obtained by taking the following kinematical condition, 

p 2 = 2p ■ q . (9) 

Note that this kinematical condition is also the consequence of the on-shell condition for the 
participating nucleons, q 2 = m 2 and (q—p) 2 = m 2 , the condition in which the physical 7rNN 
coupling should be defined. With this condition, the residues of the double poles are the 
same and in principle the double pole can provide the results independent of the coupling 
schemes. A nice aspect of this kinematical condition is that the scheme- dependent Dirac 
structure, namely 75^, which appears only in the pseudovector scheme, disappears and there 
are only three distinct Dirac structures common to both schemes. 

Now we consider 75 p 1 structure as studied by Birse and Krippa (BK). Using the pseu- 
doscalar interaction Lagrangion, the phenomenological side of this structure takes the form 

ig vN X%m + _ _ _ 



(q 2 — m 2 + ie) [(q — p) 2 — m 2 + ie] 

BK carried out the calculation using this phenomenological form with taking p^ = in the 
denominator and their sum rule formula can be succinctly written as 

g nN + AM 2 = f(M) . (11) 

f(M) is a function of the Borel mass M. The simple pole contribution associated with the 
transition N — > N* is contained in the unknown constant A. To estimate the contribution 
from A, BK first neglected A and determined g n ^. In their second method, they applied the 
operator 1 — M 2 d/dM 2 to eliminate the uncertainty associated with A. This second method 
is equivalent || to constructing the sum rule for (q 2 — m 2 )U(q,p), which was proposed by 
Jin . Anyway, the two prescriptions of BK yield the results differed by less than 5 % and 
therefore the simple pole contribution from the continuum is claimed to be negligible. 

However, the statement about smallness of A is actually scheme-dependent. If the pseu- 
dovector interaction Lagrangian is used, the phenomenological side for the 75 p structure 
takes the form 



4 



ig^N^N i i9^N^N m i n 2 \ 

2m (q — p) 2 — m? + ie (q 2 — m? + ie)[(q — p) 2 — m? + ie] 

Note that the double pole is the same as above but we have an additional simple pole of 
the nucleon which in principle could change the results. If this were used in BK's sum rule, 
then instead of Eq. flTT|), BK would have obtained 



9l Jl-- +M 2 = /(M). (13) 




Since M ~ m ~ 1 GeV, the value of the new term is about 0.5. If we neglect A as BK did 
in their first method, then this sum rule could give the result twice of what BK obtained. 
On the other hand, the second method of BK leads to the same result of BK, indicating 
that the contribution from A could be large. 

One more Dirac structure, 750"^, contains the double pole only, independent of the 
coupling scheme. The common phenomenological side for this structure is given by 



(q 2 — m 2 + ie) [(q — p) 2 — m 2 + ie] 

This structure is zero in the soft-pion limit, but beyond the soft-pion limit as BK did, this 
can provide additional sum rule for independent from the BK sum rule. Since there is 
no simple nucleon pole in this case, the simple pole structure comes only from N — > N*. 
With this in mind, we will construct the sum rule for this structure. As we will see, this 
structure has also less uncertainty in the OPE. 

Now, we construct QCD sum rule for g n N by considering •y5<7^ lJ p fJi q u structure. Specifi- 
cally, we consider the correlation function 



U(q,p) = 1 J rf 4 xe^'(0|T[J p (a;) J n (0)]|7r + (j9)) . (15) 

Here J p is the proton interpolating field of Ioffe 0, 

J P = e a bc[ulC'j ll u b }'y 5 'j' M d c (16) 

and the neutron interpolating field J n is obtained by replacing (u, d) — > (d, u). In the OPE, 
we will only keep the diquark component of the pion wave function and use the vacuum 
saturation hypothesis to factor out higher dimensional operators in terms of the pion wave 
function and the vacuum expectation value. This will be more or less similar to BK, with 
some distinctions as we will see. 

The calculation of the correlator, Eq. flT5|), in the coordinate space contains the following 
quark- ant iquark component of the pion wave function, 

D a J = (0|<(*)4(0)k + (p)> • (17) 

Here, a and (3 are Dirac indices, a and a' are color indices. The other quarks are contracted 
to form quark propagators. This matrix element can be written in terms of three Dirac 
structures, 
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£>3 = % : (7 At 75) a/3 (0|rf(0) 7 , 75 n(x)|7r + (p)) + ^(z7 5 ) a/3 (0|J(0)z 75 n(a;)|7r + (p)) 

!U -(75^) a/3 (0|J(0) 7 5^«(x)|7r + (p)) . (18) 



24 



These matrix elements can be written in terms of pion wave functions [flOf] . Since we are 



doing the calculation up to the first order of p^q u -, we need only the overall normalization 
of the wave functions. In fact, to leading order in the pion momentum, the first and third 
matrix elements are given as [[TTJ, 

(0|J(0)7 M 75w(x)|vr + (p)) = iV2U Pli + 0(x 2 ) (19) 

(0|J(0) 7 5^(x)|7T + (p)) =iV2(p (t x v - Pv x (t ) ftf S * m * , . (20) 

6(m u + m d ) 



In Eq. (|T9"|), 0(x 2 ) contains the twist 4 term whose contribution was denoted by 5 2 in 



Ref. [f|||. This term does not contribute to our sum rule up to the dimension we consider 
below. Note that in Eq. ( p0[ ) the factor f. K m 2 n; /{m u + can be written as — (qq)/f n 
by making use of the Gell-Mann— Oakes— Renner relation. Although the operator seems 
gauge dependent, it is understood that the fixed point gauge is used throughout. It is then 
interesting to note that the left hand side of Eq. ( p0|) can also be expanded in x such that 
the matrix element that contributes is effectively one with higher dimension, 

(0|J(0)75^ j D qM (0)|tt + (p)) = iy/2^g QV - Pv g^) ^J^— . (21) 

The first term in Eq. ([Tj|) is the one considered in BK's sum rule for the 75^ structure. 
However, the other two matrix elements can contribute to the 75^ structure as the calculation 
is being done beyond the soft-pion limit. Indeed, up to the highest dimension 7 that BK 
considered, the third matrix element contributes to the 75^ structure. Note that the chirality 
of the third term is different from the first term. Because of the chirality difference, the chiral 
odd condensate contributes to 75^ when the third matrix element is taken for D^,. That is, 
if we include the third term of Eq. flT8|), because it is effectively dimension 4 [ Eq. (51)], one 



dimension higher than Eq. (|T9"D Q we could have an additional odd dimensional operator, 
(qq), to form a dimension 7 operator. It only changes the highest dimensional operator in 
the BK's sum rule [Eq.(12) of Ref. Q] but its contribution is much larger than their highest 
operator because it contains the quark operator. 

In our sum rule for the 75<7 MV p At g I/ structure, the second matrix element in Eq. ([Tj|) does 
not contribute up to dimension 7. It is straightforward to calculate the OPE and up to 
dimension 7 and we obtain 



V2(qq) 



\n(-q 2 ) + ^fn_ /<^g2\ 1 + m2 ofn 



3 q 2 \ 7T / 216/ w g 4 Qq 



(22) 



1 The physical dimension of Eq. ( |i~9D is actually 2. Ordinarily, such operator in vacuum is of 
dimension 3, which is the way that BK counted the dimension of Eq. (|i~9|). In this way of counting, 
Eq. ( |21~|) is of dimension 4. 
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In obtaining the first and third terms, we have used Eq. fl20|) while the second is obtained by 
taking the first term in Eq. fll8l) for the matrix element D^, and replacing one propagator 
with the quark condensate. The fourth term is also obtained by taking the first term in 
Eq. fll8f) but in this case other quarks are used to form the dimension five mixed condensate, 
(qg s a ■ Qq), which is usually parameterized in terms of the quark condensate, m\(qq) . The 



From Ref. |TTJ, we take the range, 0.6 < m,Q < 1.4 GeV 2 



value of m 2 , is not well-known, 
and see the sensitivity of our results. One interesting aspect of our OPE expression is that 
the quark condensate, which could be the one important source of uncertainty in the final 
result, is just an overall factor. 

By matching the OPE expression with its corresponding phenomenological side after 
taking Borel transformation, we obtain 



\ 2 

9-kN^n 



B + w- 



-m? /M 2 



(gg) 



M 2 E, 



127T 2 



1 <fl 
216M 2 6M 2 



• (23) 



Here x n = S^n/M 2 with SVtv being the continuum threshold, and E n (x) — 1 — (1 +x + • ■ 
■ + x n /n\) e~~ x . The unknown parameter B comes from the simple pole of TV — > N*. To 
eliminate the parameter in the LHS, we use the chiral-odd nucleon sum rule which is 
given by 



m\ 2 N e- m2 ' M2 



M 4 E 1 (x N ), 1 



7T 



(24) 



where xn = Sn/M 2 with Sn being the continuum threshold for the nucleon sum rule. As 
BK did, we take the ratio of Eqs. (|23| ) and (^4j) to obtain 



9ttn(BM 2 + 1 




a + bM 2 

+ 16vr 2 / 2 M 2 /3 + vr 2 Q 2 ) /54 - 27r 2 m 2 / 2 /3 



M 4 E 1 (x N ) 



^ 2 )/6 



(25) 



It is now clear why our sum rule has advantages over BK's sum rule. In BK's sum rule, there 
are two main sources for the error. First one is the uncertainty in the twist-4 term, 5 2 in their 
notation, which yields ±2 error in their final result for g,^- The additional ±2 error comes 
from the uncertainty in the value of the quark condensate. These large errors are due to the 
fact that these two terms are the main terms in their OPE. No such uncertainties enter in 
our sum rule. Instead, we have different sources for the error, m 2 , and the gluon condensate. 
Note that the two sources are the highest dimensional term in our OPE. Therefore, they 
should be suppressed in the Borel window chosen. However, as m 2 , is very uncertain, the 
result has some dependence on this parameter. The error due to will be investigated 
in our numerical calculation. The additional error due to the gluon condensate, which we 
take l^-Q 2 ^j = (0.33 GeV) 4 in our numerical calculation, is very small because firstly its 
contribution to the sum rule is very small and secondly its variation is canceled in the ratio, 
Eq. 

To extract the value of g n N, we take different approach from BK. Normally, the unknown 
constant B is eliminated by applying the operator, 1 — M 2 d/dM 2 , to both sides of the sum 
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rule. It is fine mathematically but, as the OPE does not represent the total strength of 
the correlator, this operation for extracting some physical parameter is rather dangerous. 
For example, when this operation is applied to Eq. (B^) after multiplying e m ' M , the Borel 
curve for m\ 2 N is quite different from the one obtained before the operation, depending on 
the Borel mass but more or less about 20 % difference. Furthermore, the curve for g n N, if 
obtained from this operation on Eq. (|25"D , is not stable with respect to M 2 . Such g n N is just 
a decreasing function of M 2 , showing no stability . Nevertheless, g n N obtained by following 
the two methods of BK around M 2 ~ 1 GeV 2 is about 9.2 or 11.8. 

In our analysis, we directly fit the RHS of Eq. (^) with the function a + bM 2 within a 
reasonable Borel window and determine the coefficients a and b. To check the reliability of 
the fitting process, we also calculate \ 2 — Y^=\[RHS{Mi) — LHS(Mi)} 2 /N. The nucleon 
continuum Sn is set to be 2.07 GeV 2 which corresponds to the mass squared of the Roper 
resonance. By restricting the continuum contribution to be around 50 % of the first term in 
the nucleon sum rule of Eq. fl2"4"|), the maximum Borel mass is found to be around M^ ax ~ 
1.24 GeV 2 . The minimum Borel mass is determined by requiring the highest power correction 
in the pion sum rule of Eq (|23| ) to be less than 10 % of the total OPE. Such chosen window 
provides the common region of the two sum rules. For example, using m 2 , = 1 GeV 2 , this 
condition leads to M^ in ~ 0.82 GeV 2 . We take SVjv to be equal to Sn in our analysis. 
But the sensitivity to S n N is very small. For S^n ~ Sn = ±0.5 GeV 2 , g n N changes only by 
±0.3. Within the Borel window, the parameters, a and b, are determined by the best fitted 
method. Then, the coupling, g n N, is obtained via g n N = mo,/ f n using the physical values, 
m = 0.94 GeV and U = 0.093 GeV. 

The results are listed in Table | for various values of m\. First note that the contribution 
from the unknown single pole, N — > N*, represented by b is relatively small. At m 2 , = 0.6 
GeV 2 , this contribution is almost negligible. Our result is g^N — 9.76 ± 2.04 which is rather 
smaller than its empirical value of 13.4. It is however possible to obtain a larger value 
by reducing m 2 ,, by relaxing the 10 % restriction for M^ in or by restricting the continuum 
further. Any of these moves the Borel window to smaller Borel masses so that the Borel curve 
contains more contribution from higher nonperturbative operators. So our result should be 
interpreted within the standard prescription of the conventional QCD sum rule. 

In summary, we have constructed QCD sum rule to determine irNN coupling constant. 
We noted that the previous calculations have the coupling scheme dependence and the 
determined value of the coupling suffers from large uncertainties in QCD parameters. We 
have discussed the issue related to the coupling scheme dependence and stressed that one 
should look into the double pole beyond the soft-pion limit in the QCD sum rule study of 
the 7rNN coupling. In this work, we have proposed to look at the Dirac structure, 750"^, 
beyond the soft-pion limit. This structure is, first of all, independent of the effective models 
employed in the phenomenological side and further provides the irNN coupling with less 
uncertainties from QCD parameters. It would also be interesting to purse similar work with 
a different nucleon interpolating field. 
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TABLES 



TABLE I. The best fitted parameters for a, b and g^N at given values of m\. For a given value 
of mQ, M^ in is chosen such a way that the highest dimensional operators contribute less than 10 
% of the total OPE. M^ aa , is chosen so that the continuum contribution is less than 50 % of the 
first term of the OPE in the nucleon sum rule. Such chosen M^ ax is 1.24 GeV 2 . Then, a and b 
are determined from the best fitting method within the window. is obtained by using physical 
values for nucleon mass and the pion decay constant, m = 0.94 GeV and = 0.093 GeV. 



ml (GeV 2 ) 


Ml m (GeV 2 ) 


a 


b (GeV" 2 ) 


x 2 


9nN = ma/ f n 


0.6 


0.55 


1.17 


0.05 


6.3 x 10" 4 


11.8 


0.8 


0.69 


1.01 


0.18 


1 x 10" 4 


10.22 


1.0 


0.82 


0.9 


0.27 


1.7 x 10~ 5 


9.14 


1.2 


0.94 


0.83 


0.32 


2 x 10~ 6 


8.35 


1.4 


1.05 


0.76 


0.36 


2 x 10~ 7 


7.72 
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